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Abstract 

We present our results for the NNLL virtual corrections to the matrix elements of the operators 
Oi and O2 for the inclusive process b — > s£~^£~ in the kinematical region > 4m^, where is the 
invariant mass squared of the lepton-pair. This is the first analytic two-loop calculation of these 
matrix elements in the high q'^ region. We give the matrix elements as an expansion in mc/mh and 
keep the full analytic dependence on g^. Making extensive use of differential equation techniques, 
we fully automatize the expanding of the Feynman integrals in nifJmb. In coincidence with an 
earlier work where the master integrals were obtained numerically 1[, we find that in the high q'^ 
region the as corrections to the matrix elements {si^i~\Oi^2\b) calculated in the present paper lead 
to a decrease of the perturbative part of the g^-spectrum by 10% — 15% relative to the NNLL result 
in which these contributions are put to zero and reduce the renormalization scale uncertainty to 
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I. INTRODUCTION 



Flavor-changing neutral currents play an important role in the indirect search for new 
physics. For inclusive decays there exists the framework of operator-product expansion, 
which makes theoretically clean predictions possible. Of special interest in this context is 
the decay mode B — ^ Xsi~^i~. In the regions where the lepton invariant mass squared 
is far away from the cc-resonances, the dilepton invariant mass spectrum and the forward- 
backward asymmetry can be precisely predicted. 

The status of the calculation of these observables is the following: The leading logarithmic 

n 

(_LL) and the next-to-leading logarithmic (NLL) QCD contributions were calculated in [2|, 
3, ^. Next-to-next-to- leading logarithmic (NNLL) corrections to the Wilson coefficients at 
the matching scale fi ~ rriw, which required to perform two-loop matching calculations of 
;he full standard model (SM) theory onto the effective theory, have been worked out in 

3, y, 0, 3 



at the low scale 
obtained in 



The anomalous dimensions matrices needed to obtain the Wilson coefficients 
mij jrequiring up to three-loop calculations for certain entries) were 



11 



12l |. NNLL QCD corrections at the level of the matrix elements 



of the operators involved were calculated for the dilep 
for the forward-backward asymmetry in [l|, 14, isl . 
order 1/m^, 1/m^ and 1/m^ have been worked out in 



;on invariant mass spectrum and 



m 
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17 [.Power corrections of the 



20 



21 
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2^. Finally, 



27l | certain classes of logarithmically enhanced electromagnetic corrections were 



taken into account. 

So far, analytic results for the NNLL QCD corrections to the matrix elements associated 
with the operators Oi and O2 are only available in the region of low o^. The corresponding 
results were obtained as a double-expansion in rric/mh and q^/m1 [isl. [l^. 



15|. The present 

paper deals with the NNLL QCD corrections in the high region, i.e. q^ > 4m^. In 
particular we evaluate virtual QCD corrections to the matrix elements of the operators Oi 
and O2 at order a^. In contrast to {l|, where the relevant master integrals were calculated 
numerically, we present these matrix elements as analytic functions of rric and q^. The 
purpose of the present paper is twofold: First, to deliver a non-trivial independent check of 
the results found in [l| and second, to provide the user with analytic formulas in which the 
parameters {mc/mh and fi/rrify) and can easily be changed. 

To get these analytic results, we perform an expansion in mc/mi, and keep the full analytic 
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dependence on o^. We expand the two-loop Feynman integrals by combining method of 



regions 
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311 ] and differential equation techniques 
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35|. We end up with 



an expansion of {si~^i~\0i^2\b) up to the 20th power in mc/mi,. As the resulting expressions 
for these matrix elements are rather lengthy, we are not able to print them in the paper. We 
provide Mathematica and C++ code of our results in the source files of the present paper at 
arXiv. 

The well-known breakdown of the A/m;, expansion at the endpoint = m\ seems to 
question the relevance of the perturbative contributions in the high g^-region calculated in 



this paper. However, as it was shown in [36|] and [37| (illustated there for the analogous 
lepton invariant mass spectrum in the inclusive semileptonic decay B XuC-v) that the 
integrated high g^-spectrum allows for a modified version of the heavy-quark expansion (the 
so-called hybrid expansion), our present work is well-motivated. 

The paper is organized as follows. Sections [III and IIIII are dedicated to the technical 
details of the calculation. We give all necessary definitions in Section [Tll In Section IIIII we 
explain the evaluation of the Feynman integrals in detail. In Section IIVI we investigate the 
(numerical) stability of the expansion in rric/mb, concluding that retaining terms up to the 
20th power in mc/mi, leads to precise results. In this Section we also discuss the numerical 
impact of our calculation on the dilepton invariant mass spectrum. In coincidence with [l| we 
find that in the high region the order corrections to the matrix elements {s£^i~\0i^2\b) 
calculated in the present paper lead to a decrease of the perturbative part of the g^-spectrum 
by 10% — 15% relative to the NNLL result in which these corrections are put to zero and 
reduce the renormalization scale uncertainty to ~ 2%. 



II. DEFINITIONS 

As in the previous paper [l^ we write the effective Hamiltonian that contributes to 
B Xsi+i- in the form 

„ _ 4G 



^ 1=1 
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b Oi 2 s 

— ^ — g> 




FIG. 1: Diagrams that have to be taken into account at order a^. The circle-crosses denote the 
possible locations where the virtual photon is emitted (see text). 

where we have neglected the small CKM combination V*gVub- The operator basis is defined 
as 



Oi = 




O2 = 




03 = 




O4 = 




05 = 




Oe = 




07 = 




Os = 




09 = 




Oio = 


fJ(^L7A)E.(W), 



where the subscripts L and R refer to left- and right- handed components of the fermion 
fields. The ingredients to obtain the Wilson coefficients Cj at the scale /i of order can be 
found e.g. in [3, [sl, Q . 

In the present publication we calculate the virtual a^-corrections to the matrix elements 
of Oi and O2 in the large region. Using equations of motion, we write these a^-corrections 
in the form^ 

(s£+r |0,|6)2.1oops = - (^)' [i^'(Or)trcc + i^'(09)trce] • (3) 



^ Note that because of the extra factors 1 / g1 in the definition of O7 and O9 ^ is indeed of order as . 
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The diagrams that contribute at order to 6 — sl^l are shown in FiglH By definition 



3, 



we 



we include in F^^'^ only the contributions from the diagrams in Fig. [T^-e. As in 
absorb the contribution from Fig. [If into a modified Wilson coefficient Cg. This procedure 
is convenient, because only the diagram Fig. [T]F contains infrared divergences. 

The ultraviolet renormalization works analogously to [l^. In particular we use the same 
evanescent operators. We use on-shell renormalization for the s- and the 6-field and renor- 
malize rric in the pole mass scheme. 

The kinematics is defined as follows: We denote the momentum of the incoming 6-quark 
by p and the momentum of the virtual photon by q. The momenta of the external fermions 
are on-shell such that = ml and {p — qY = 0, because we neglect the strange-quark mass. 
Furthermore we use the notations 

2 2 

s = — ^ and z = — ^. (4) 

m^ ml 

III. CALCULATION OF THE MASTER INTEGRALS 

In the present section we explain for every diagram appearing in Fig. [T] the way we 
evaluated the master integrals that are specific to it. In Appendix |X] we list the master 
integrals that appear in more than one diagram and which are straightforward to calculate. 

We use the following notation 

where = 4 — 2e. 

For simplicity we set m^ = 1 in the calculation of the master integrals, such that = z. 
The dependence of the master integrals on m?, can be easily restored by dimensional analysis. 



A. General remarks about calculation techniques 



The Feynman integrals appearing in our calculation have been reduced to a set of master 
integrals using the following methods: Tensor integrals i.e. integrals containing Lorentz 
indices have been reduced to scalar integrals via the Passarino-Veltman reduction scheme 
381]. Finally these scalar integrals can be further reduced by integration by parts (IBP) 



identities [3^ 



In particular we used the algorithm described in 



411 . To this end we 



used the Maple implementation AIR [4^ and a Mathematica implementation developed by 
us. Since we consider the region s > iz, we expanded the master integrals in z and kept the 
full analytic dependence in s. 
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For pow er expanding Feynman integrals we use a combination of method of regions 28 



30 



3l| and differential equation techniques [32, l33|, l34, l35|, |43|]. We consider a set of 
Feynman integrals Ji, . . . , /„ that depend on the expansion parameter z and that are related 
by a system of differential equations: 

= KpIfS + ga- (6) 

We obtain ([6]) by differentiating la with respect to z and applying IBP identities, from 
where we obtain the original set of integrals and further integrals contained in ga, which are 
simpler than /„ and have been calculated before. Expanding the objects appearing in ([6]) 
in e, z and In z 



la = J2 li^fe^z^ (In z) 



,0') A^J 

al3,V 

9a = (7) 

and inserting ([71) into ([6]) we obtain algebraic equations for the coefficients I^'P 

o = u + + 1)^^'"^^ - E E E ^S'^i"^ - al^f- (8) 

13 i' j' 

By means of ([H]) we can reduce higher powers in z of Jq, to lower powers. In practice this 
means that we need the leading power and sometimes also the next-to-leading power of 
as initial condition for ([H]). We have calculated these initial conditions by method of regions. 
Every region except the hard region leads to logarithms in z. As we obtain the logarithms 
occurring at leading power both from method of regions and from the recurrence relation 
([H]), differential equations provide a non trivial check for method of regions, i.e. we can make 
sure not to have forgotten or counted twice any region. 

In ([7]) we did not specify which values the summation index j takes. Indeed we will have 
to deal with integrals that come with half-integer values of j i.e. they have to be expanded 
in ^/z. On the other hand we have to presume that there exists /cmax such that laf' = for 

6 



all k > /cinax in order to solve ([H])- We use the algorithm that was described in [3J] to get the 
possible values for j and to determine fcmax- In addition this algorithm allows us to evaluate 
the coefficients ll^'f' numerically. We used this feature to test the initial conditions. 

In the following we will show in detail how to evaluate the master integrals occurring 
from the diagrams in Fig. [T] by this procedure. 



B. Diagrams of Fig. [1^ 



The topology of Fig. [Bi contains in addition to (jM]), (lAlip . (1A12I) and (1A13I) . which are 
easy to evaluate, these two master integrals 

1 



hi = j [dk][dl] 



{k + p- qy{k + pY{{k + If - z){P - z) 



where we use the notation and assume implicitly that every denominator contains a 
positive imaginary part +iO. We need both integrals in leading power i.e. at There 
are three regions that contribute to this power: The hard region k'^,1'^ ~ 1, the soft region 
/c'^ ~ 1, and ~ ^/z and the collinear region where both k and / are collinear to p — q 
(scaling see below). Both integrals get a leading power contribution in the hard region. The 
hard region corresponds to setting 2; = in the integrand. In this limit we can reduce Ia2 
to lai by IBP identities, hi at z = can by evaluated via Feynman parameterization to 

/ - ^ . 2^..+..^2. ^(6)^(26)^3(l-e)^(l-26) 

) r(i + 6)r(2 - 26)r(2 - 3e) '^'^^'^ i' i + ^. i - m 



where 



,F,{a, b; c; x) = dtt'-\l - tr~'~\l - tx)"^ (11) 



with 3?c > > 0. We used the Mathematica packages described in [4^, |45| to obtain the 
expansion in e of 2F1. 

In the soft region /c^ ~ 1, ~ ^/z only Ia2 gets a leading power contribution: 

= I {k+p-qmlpwip-ir 

Using IBP identities, (fT2l) can be reduced to a product of two simple one-loop integrals. 
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Let us consider the collinear region. We introduce the following light-like vectors and 



n_, which fulfil 



n 



and n+ ■ n_ = 1. We define the decomposition of a Lorentz 



vector into light-cone coordinates: 



(13) 



where = k-n^:. We choose n+ to be collinear to p — q and introduce the following scaling 



k±,l±r^^/z and k ,1 ^ z. 

As before, only Ia2 gets a leading power contribution in this region. 

1 



'a2,c — Z 



-2e 



[dk][dl] 



{k+p- qy{2k+p~ + 1)((A; + If - 1){P - if 
Via Feynman parameterization we evaluate ( fT5l) . obtaining 

-1 



2e r2(e) 



2F1 l,l;2-e;l-s . 



(4vr)4 V- " ^ 2r(l - e) 
Finally the leading power contributions of the master integrals up to order read 



(14) 



(15) 



(16) 



r(0) 
-'al 



.4e 



1 

2^ 



In (s) — -ivr — I 



(47r)4' 

In^ (s) + (5 + 2i7r) In (s) + Lis (1 - s) 



12 2 



r(0) 
'a2 



(4vr) 



In (s) In^(^) + (ivr In (s) + Lia (1 - s)) \n{z) 



In (s) + ivrLis (1 — s) — Lis (1 ~ s) 



(17) 



We continue with the calculation of the sub leading powers of lai and Ia2- By differenti- 
ating lai and Ia2 with respect to z and applying IBP identities we obtain a coupled system 
of differential equations of the form ^ with h starting at order e*^ and z~^. More explicitly 
([H]) becomes 

«+2 j 



o = (j + i)/lS'''^ + (fc + i)/, 



\^ \^ lO") rO-i',fc) _ 

/ J / J / J "'af3,i''^al3,i-i' ya,i ' 
(3=1,2 i'=0 j'=-l 



(18) 



From ( fTSl) together with (fTTI) we obtain the sub leading powers in z of lai and /qs- We 
also obtain the coefficient of the z^lnz-term of Ia2, which we already calculated in (fTTI) . 
This means that the differential equations provide a non-trivial check for method of regions, 
which was used for the leading power calculation. 
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C. Diagrams of Fig. [TJ) 



The topology Fig. [TJd comes with the master integrals 

1 



= J [dk][dl] 
h2 = J [dk][dl] 



{{k + py- 


mk + p- 


qf- 

1 


l){P-z){{k^ 


-ly-z) 


{{k + py- 


mk+p- 


qf- 

1 


l)(/2-^)2((A;- 





hz J [dk][dl] ^ _ + p _ _ i)2(p _ ^)((^ + ly _ ^) • (19) 

We need these integrals in leading power. Besides the hard region, where all of these integrals 
get a leading power contribution, Ii,2 also gets contributions from two further regions. In 
the soft region defined by /c'^ ~ 1 and l'^ ~ y/z Ib2 becomes 

'"'^ = ""^ /[^^] [^^1 {{k + py-im + p-qy-iW-iyk^^ 

which is a product of ( ]A8I) and a trivial tadpole integral. In the coUinear region defined by 
/c'^ ~ 1, Z"*" ~ 1, Z_L ~ ^/z and /~ ~ z, 1^2 takes the form 



However, the collinear region has an overlap with the soft region, where fl2ip reduces to 
On the other hand fl21l) is indeed equal to fl^Ul) which can be seen by the following argument: 
Consider the integration [dl]. The integrand depends besides on terms constant in only 
on P and = I - n^. So n't is the only Lorentz vector that multiplies l'^. Because of Lorentz 
invariance the integral can only depend on through = 0. So we can set to zero 
such that (12T!) reduces to (l20l) . This is to say the collinear region has already been taken 
into account by the soft region. To avoid double counting we have to skip the contribution 
([21]). Analogously we can introduce another collinear region /c^ ~ 1,/ ~ By the same 
argument we see that also this region has been already taken into account in ( !20l) . 
In the hard region IBP identities provide a reduction of ( |T9l) to the set of integrals 

1 



hi,h = j[dk][dl] 



{{k + pY -l){{k+p- qf - l)P{k + /) 



2 



h2,K - / [dk][dl] ((^ + p)2_l)((^+^_g)2_l)^4(^ + ;)2- (22) 



These integrals can be evaluated via differential equations with respect to s. By defining 

hi,h 

h2,h 



(23) 



and differentiating / with respect to s we obtain a differential equation of the form 

d 



ds 



I = hl + g 



(24) 



where g contains the integrals (1A14I) and (1A15I) . We define the expansion in e 

oo 

i=-2 

oo 

h = 5^/i«e^ 



i=Q 

oo 



9 = 



(25) 



i=-2 



and write ( 1241) in the expanded form 

d 



j(-2) ^ ^(0)/(-2) ^ ^-2) 

(is 
as 
ds 



(26) 



In our special case h^i2 — such that ( !26l) decouples and we can solve (|26l) by the common 
methods separation of variables and variation of constants. From Feynman parameterization 
we see that the limit s = does not lead to additional divergences in e and can be used as 
initial condition for 0261) . Finally we obtain 



hl,h 



(4vr) 



4e 



-1 2J^arcsin^-f 3 (s - 3) arcsin^ ^ 

\ 1 — 

2e^ e s - 1 



-Stt^s - 114s + Ttt^ + 114 4:- SI ys , , , 
12 (s- I) + V ^ 1~ (4 - s) + In (s) + 10) 



+ CI2 I 2 arcsin j — 2CI2 ( 2 arcsin ^ — h tt 



.4e 



(4vr) 



6 arcsin 



2 7r2 



^12 In (1 - s) + 3 In (s)) arcsin^ — 



V S \ V 5 / V 5 

4CI2 I 6 arcsin h ir arcsin In (1 — s) — 301^^ 2 arcsin 

' 2 / 2 3 ^ ^ I 2 



+ 6CI3 ( 2 arcsin ^ + ] + ^Clg ( 6 arcsin ^ + vr ) + 3C(3) 



(27) 



where Cl2(0) = ^5Li2(e*'^) and CW) = ^Li^^e 



i4>\ 
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D. Diagrams of Fig. \Vp 



The topology Fig. [T]3 comes with the master integrals 

1 



Id 

/c3 



[dk][dl] 

[dk][dl] 
[dk][dl] 



{P- 


z){{k + ly - z){{i ^ 
1 




- z){k + p- 


■qY 




zf{{k^lY-z){{l^ 

1 




- z){k + p 


-qy 



(28) 



(/2 - z){{k + If - zY{{l + qy - z){k + p- qf ■ 
They all get leading power contributions from the hard region, where IBP identities lead to 
a further reduction of Ic2 and Jcs to Id- Id can be calculated by a differential equation with 
respect to s, which reads: 

1 



d. 2s- 1 ^ 



where 



fi:(/i,e) 



(47r)4 ' 



2i-Ke 



T{-l + 2t)T\l-e) 
r(2 - 3e) 



The most general solution of (!29|) is given by 
cs-%1 - s)-^ - K(/i,e)(l - sy 



2Fi(-e,e;l + e;s) ^ g2Fi(l - e, 1 + e; 2 + e; s) 



(29) 



(30) 



(31) 



e 1 + e 

where we have to determine c. We note that both s = and s = 1 are no appropriate 
initial conditions. Hence we determine c by calculating the term proportional to in the 
expansion of Ici,h around s = 0. The Mellin-Barnes representation (see e.g. 3l|]) of Id^h 
reads 



'cl,/i 



(/i^e 



>2e 



(4vr) 



/ioo pi 
dts^T{~t)T{t + 2e) dxx~'^'{l-x) 
■ioa Jo 



X [' d'yy^'-^-\l - y,)-^y^'^-\l - y^y^l - ym) 
Jo 



(32) 



We have to calculate the residue at t = — e in (132|) . which arises due to the integration 



1 ^2„,„,-l-e-t/ 



lo d^yy'i 



at yi = 0. So we can set ?/i = in the ellipsis and obtain 



IclM = s' 



(47r)4 (1 - 2e)r(2 - 3e) 

where the ellipsis denotes integer powers of s. Hence c reads 



(33) 



(/i^e 



,2€ 



(47r) 



;i - 2e)r(2 - 3e) 



(34) 
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In the collinear region /c"*", /+ ~ 1, k±,l± ~ y/z, k ,1 z both 1^2 and Jc3 get a leading 
power contribution: 

1 



Ic2,c = z-^' j [dk][dl] 



(/2 - 1)2((A; + /)2 _ l)(2/+g- + + 2A;-(p - 



(47r)2 2(1 - e) s ' ' V ' ' ^ 

1 



h3,c = z-'' l[dk][dl] 



(/2 - 1)((A; + /)2 - l)2(2/+g- + s){k^ + 2A;-(p - 



(47r)2 2(1 - e) s " ' V ' ' ' ^ 

where 3F2 is given by 

3F2(ai,a2,a3;6i,62;a;) = > p. s 7^77— — xp/, , — VH"' 

^ T{ai)T{a2)T{a3) T [bi + n)T {b2 + n) nl 



and can be expanded in e by the tools developed in 4^, |45| 



In the soft region fc'^ + ~ ^/z only Jc3 contributes in leading power: 

= ""^ /[^^] /2(/ + g)2(/-p + ,)2(p_,)2' (37) 
which is a product of two simple one-loop integrals. There are two further collinear regions 
A; + Z ~ 77,+ and k + I ~ n^, where Jcs obtains a leading power contribution. However 
by an argument similar to that given in the previous subsection we can show that these 
contributions have already been taken into account by flHTj) . 

As described above the subleading powers of ( |28l) are obtained via differential equations 
with respect to z. Like in the previous cases the terms of the order In z provide a check 
that we have taken all regions contributing at leading power consistently into account. 

E. Diagrams of Fig. [T]i 

The topology in Fig. [T]i comes with two sets of master integrals. 



-^dll = 


j [dk][dl] 








1 








{P- 


Z){{1 


-kf- 


m- 


q?- 


z){{k- 


-p)2-l) 


^(112 = 


1 [dk][dl] 








1 








(/2- 




-ky- 


zm- 




-z){{k 




Idl3 = 


1 [dk][dl] 








1 








{P- 




-kf- 


m- 


g)2- 


2)2((A; 


-p)2-l) 


Idli = 


1 [dk][dl] 








1 








(F- 




-ky- 


z){{l- 


qf- 


z){{k~ 


-p)2-l)2 



(38) 
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and 



Id2i = j[dk][dl] 
h22 = j [dk][dl] 



k^{{l- 


kf- 




-q?- 
1 


z){{k~pY-i) 


k^{{l- 


ky- 




-qf- 
1 


zY{{k-pY-i) 



/.23 - j [dk][dl] p((;_^)2_,)((^_^)2_,)((fc_p)2_l)2- (39) 

Let us consider the first set fl38p . In the hard region this set reduces by IBP identities 
to Jrfii and Idi2- These integrals can be calculated by differential equations with respect to 
s. We obtain a system of differential equations similar to ( 1261) where we have to use s = 1 
as initial condition because the integrals diverge at s = 0. The matrix h^^^ has vanishing 
off-diagonal elements such that the system of differential equations decouples. In addition 
the /i^*-* contain only terms of the form 1/(1 — s), 1/s and s". So we can reduce the integrals 



to harmonic polylogarithms, which were defined in 



461 ] ■ The way to do this is very well 



491] to convert 



described in Section 2.4 of 47|]. Finally we used the program described in 48 
harmonic polylogarithms into common functions like polylogarithms. 

The soft region — k^ ^/z leads to a leading power contribution of 1^12 

'''''' = I k-{k-qmk-py-iw-ir 

where we substituted I ^ I + k. This integral is a product of (1A5I) and a simple one- loop 
tadpole integral. 

The soft region — ~ ^/z leads to a leading power contribution of Idiz 

z~' r „ 1 



(41) 



s J ' ' A;2(P-l)(/2-l)2' 

where we substituted I ^ I + q. This integral is a product of two simple one-loop integrals. 

Let us consider the second set of master integrals (|39l) . In the hard region the set reduces 
via IBP identities to 1^21 • We evaluated Id2i by a differential equation with respect to s. 
Solving this differential equation is a straightforward calculation, which is analogous to the 
way we solved (l29l) . 

The soft region — q^ ^ ^/z leads to a leading power contribution of 7^22 

= /[^^] k^ik-qyuk-pr-iw-ir ^''^ 

which coincides with f HUj) . 
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Besides the leading power we also need the order z of ld2\- It is straightforward to 
calculate the order z contribution of the hard region by expanding the integrand of ld2\ up 
to the order z. Finally the soft regions — k^ ^ -^i and — ^ ^ contribute at order 
z. Since these regions do not overlap we have to take both of them into account. After an 
appropriate shift of /, ld2\ can in both regions be cast into the form 

= I k^ik-qmk-p)^-iw-iy 

which is similar to (HOl) . 



F. Diagrams of Fig. [T^ and f 



The integrals occurring in the diagrams of Fig. [T^ reduce to ( ]A2p . ( lAllI) . ( ]A12[) and 
f lAlSp . The topology of Fig. [Tf factorizes trivially into two one-loop topologies, which have 
already been evaluated exactly in s in IJ]. As already mentioned in Section [Tll Fig. [Tf does 



not contribute to the form factors by definition; its effect is, however, absorbed into a 



modified Wilson coefficient Cg as in 14j ]. 



IV. RESULTS 

A. Results for the form factors F^^2^^ in the high region 

In Section IIIII we calculated the two-loop diagrams in Fig. [1^ -e which contribute to the 
form factors -Fi 2'^"* defined in dH]). In addition, there are counterterm contributions which have 
to be taken into account. These counterterms are qualitatively the same as those discussed 
in Section III.B of [l^. Because its calculation in the high g^-region is straightforward, we 
do not list their explicit results. We only stress that in the following results the c-quark 
mass is renormalized in the pole-scheme. 

We calculated the renormalized form factors -Fi 2'^"* in the large g^-region as expansions 
of the form Cnm{s) z"' lia^ z {n = 0, m = 0,1,2,...), keeping the full analytic 

dependence on s {z = ml/ml, s = /ml). We included all orders up to 2^°. To demonstrate 
the convergence of the power expansions, we show in Fig. [2] the form factors as functions of 
s, where we include all orders up to z^, z^ and 2^°. We use as default value z = Q.l such 
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FIG. 2: Real and imaginary parts of the form factors functions of s. To demonstrate the 

convergence of the expansion in z we included all orders up to z^, z^ and z^^ in the dotted, dashed 

and solid lines respectively. We put /i = ruh and used the default value z = 0.1. 
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s 






0.25 


0.6 


-0.928 - 0.408i - 0.856^ 


5.57 + 2.45i + 5.14£ 


0.7 


-0.909 - 0.458i - 0.856^ 


5.45 + 2.75i + 5.14£ 


0.8 


-0.888 - 0.500i - 0.856^ 


5.33 + 3.00i + 5.14£ 


0.9 


-0.867 - 0.535i - 0.856£ 


5.20 + 3.21i + 5.14£ 


0.27 


0.6 


-0.919 - 0.347i - 0.856£ 


5.52 + 2.08i + 5.14£ 


0.7 


-0.905 - 0.402i - 0.856£ 


5.43 + 2.41i + b.Ue 


0.8 


-0.888 - 0.449Z - 0.856^ 


5.33 + 2.69i + 5.14^ 


0.9 


-0.869 - 0.488i - 0.856^ 


5.21 + 2.93^ + 5.14^ 


0.29 


0.6 


-0.904 - 0.280i - 0.856^ 


5.42 + 1.68i + 5.14^ 


0.7 


-0.896 - 0.342i - 0.856^ 


5.38 + 2.05i + 5.Ue 


0.8 


-0.883 - 0.393i - 0.856£ 


5.30 + 2.36i + 5.14£ 


0.9 


-0.867 - 0.437i - 0.856£ 


5.20 + 2.62^ + 5.14^ 


0.31 


0.6 


-0.879 - 0.208i - 0.856£ 


5.28 + 1.25i + 5.14£ 


0.7 


-0.881 - 0.277i - 0.856^ 


5.29 + 1.66i + 5.14£ 


0.8 


-0.874 - 0.334Z - 0.856^ 


5.24 + 2.00i + 5.14^ 


0.9 


-0.862 - 0.382i - 0.856^ 


5.17 + 2.29^ + 5.14^ 


0.33 


0.6 


-0.842 - 0.130i - 0.856^ 


5.05 + 0.779i + 5.14^ 


0.7 


-0.858 - 0.207i - 0.856£ 


5.15 + 1.24^ + 5.14^ 


0.8 


-0.859 - 0.269i - 0.856£ 


5.15 + 1.62f + 5.14^ 


0.9 


-0.853 - 0.322i - 0.856£ 


5.12 + 1.93i + 5.14£ 



TABLE I: Numerical results for the form factors i^x 2 , 



for different values of z and s (i = ln—' 

^ nib ' 



that the cc-threshold is located at s = 0.4. One sees from the figures that far away from the 
cc-threshold, i.e. for s > 0.6, the expansions for all form factors are well behaved. 

In Tab. [Ifni we list numerical values of the form factors for different values of z and 
s, retaining the dependence on the renormalization scale fi. We compared our values in 



Tab. [mm with the numerical values 



50l ] that were used in [l|. We obtain nearly perfect 
agreement, i.e. the difference is always smaller than 1%. 

Unfortunately, the form factors are too lengthy to be given explicitly in this paper. 
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s 






0.25 


0.6 


8.72 - 22.9i + (-5.47 - 3.23i)^ - 1.05^^ 


13.2 + 13.6* + (22.1 + 19.4«)£ + 6.32^^ 


0.7 


9.53 - 19.8i + (-5.13 - 3.31i)£ - 1.05^^ 


10.3 + 14.5« + (20.1 + 19M)1 + 6.32^2 


0.8 


9.92 - 17.5i + (-4.86 - 3.36z)^ - 1.05^^ 


7.88 + 14.8i + (18.5 + 20.2i)£ + 6.32^^ 


0.9 


10.1 - 15. 8i + (-4.64 - 3.40i)£ - 1.05^^ 


5.93 + 14.8z + (17.2 + 2QAi)i + 6.32f 


0.27 


0.6 


7.65 - 26.6i + (-5.66 - ?>.\\i)i - 1.05^^ 


15.0 + 10.9i + (23.3 + 18.7i)f + 6.32^^ 


0.7 


9.07 - 22. 7i + (-5.29 - 3.23i)^ - 1.05^^ 


12.0 + 12. 6i + (21.1 + 19Ai)e + 6.32^^ 


0.8 


9.78 - 20.0i + (-5.00 - 3.30i)^ - 1.05^^ 


9.44 + 13.3« + (19.3 + 19.8i)£ + 6.32f 


0.9 


10.2 - 17.9i + (-4.76 - 3.35i)^ - 1.05^^ 


7.35 + 13.6« + (17.9 + 20.1i)£ + 6.32^^ 


0.29 


0.6 


5.76 - 31.0i + (-5.88 - 2.95z)^ - 1.05^^ 


16.6 + 7.46i + (24.6 + 17.7zy + 6.32^2 


0.7 


8.11 - 26.2i + (-5.47 - 3.12z)£ - 1.05^^ 


13.6 + lO.H + (22.2 + 18.7i)i + 6.32^^ 


0.8 


9.32 - 22. 8i + (-5.15 - 3.22i)£ - 1.05^^ 


11.0 + 11. 5i + (20.3 + 19.3i)e + 6.32^^ 


0.9 


9.98 - 20.3i + (-4.89 - 3.29i)£ - 1.05^^ 


8.81 + 12.2i + (18.7 + 19.7i)e + 6.32f 


0.31 


0.6 


2.65 - 35.9i + (-6.12 - 2.74i)£ - 1.05^^ 


17.9 + 3.06i + (26.1 + 16Ai)i + 6.32^2 


0.7 


6.46 - 30.1i + (-5.67 - 2.98i)£ - 1.05^^ 


15.1 + 7.05i + (23.4 + 17.9i)i + 6.32^^ 


0.8 


8.41 - 26.0i + (-5.32 - 3.12i)^ - 1.05^^ 


12.5 + 9.24i + (21.3 + 18.7i)£ + 6.32P 


0.9 


9.50 - 23.0i + (-5.04 - 3.21z)^ - 1.05^^ 


10.3 + 10.5« + (19.6 + 19.3i)£ + 6.32^^ 


0.33 


0.6 


-2.28 - 41.7i + (-6.39 - 2.45i)^ - 1.05^^ 


18.4 - 2.6H + (27.7 + U.7i)i + 6.321'^ 


0.7 


3.84 - 34.6i + (-5.89 - 2.79i)£ - 1.05^^ 


16.3 + 3.20i + (24.7 + 16.8i)i + 6.32f 


0.8 


6.90 - 29.7i + (-5.51 - 2.99i)£ - 1.05^^ 


13.9 + 6.42i + (22.4 + 17.9i)i + 6.32f 


0.9 


8.60 - 26.1i + (-5.20 - 3.12i)^ - 1.05^^ 


11.7 + 8.3H + (20.5 + 18.7i)e + 6.32^^ 



TABLE II: Numerical results for the form factors 2 , foi^ different values of z and s (^ = In ^). 

Hence, the complete analytical results are attached to the source-code files of the present 
paper at www.arxiv.org. The Mathematica file FJiigli.m contains the expressions for 
FlTHighRe, F17HighIni, F19HigliRe, F19HighIni, F27HighRe, F27HighIm, F29HighRe and 
F29HighIm, which represent the real and imaginary part of the form factors -F^ 2'^'' in the 
high s region; they are defined in terms of muh, z and sh standing for fi/mf,, z and s respec- 
tively. Additionally this file contains the expressions for DeltaF19HighRe, DeltaF19HighIin, 
DeltaF29HighRe and DeltaF29HigliIin, which have to be added to the pole-scheme form fac- 
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' ' ' ' ' ' ' 

0.6 0.65 0.7 0.75 0.8 0.85 0.9 

s 

FIG. 3: Perturbative part of R{s) at NNLL. The solid line represents the NNLL result, whereas in 
the dotted line the order as corrections to the matrix elements associated with Oi^2 sue switched 
off. We use n = rrij,. See text for details. 

tors in order to switch from the pole-scheme to the MS-scheme of the c-quark mass. For 
completeness we also provide the file F_low.m, which contains the analogous expressions in 
the low s region (FlTLowRe etc.) taken from [l^. For numerical purposes we also provide 
the C++ header files F_l .h and F_2.h that contain the analogously defined functions 

double F_17re (double muh, double z, double sh) , 
double F_17im(double muh, double z, double sh) , etc. 

valid in both high and low s region. These files need for numerical evaluation of the harmonic 
polylogarithms the header file hpl.h, which we provide at the same place. 



B. Impact on the dilepton invariant mass spectrum in the high region 

In this section we briefly discuss the impact of the form factors -Fi 2'^'* calculated in this 
paper on the g^-spectrum at high values of q^. To this end, we consider as in jl^ the 
perturbative part of the ratio 

Ms) = I ^r(i^ ^ xjn-) 

^^'^ r{B^X,e-u,) ds ' ^^^^ 
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where the formulas for the decay rates r{b X^e v^) and dT{h Xsi~^£ ) / ds can be found 
e.g. in Section VI of jl4|. The parameterization of dT{b Xsi~^i~) /ds as specified in (89) 



and (90) of [l^ is also valid in the high region. All the ingredients contained in these two 
eqs. are available for arbitrary g^, except -Fi 2 8'- The expressions for -F^^2^'' were derived in 
the previous sections of this paper in the high range. The calculations of the renormalized 
form factors Fg^'^'' is much easier and we therefore immediately give the results (valid for 
arbitrary s G [0, 1]): 

~ 27(1-5)4 9 (1-S)2 9(1-5)3^^ bs + Zs )^Tcsmy^ 

16 2 + s . 2(^\ 8s 32 u 8 . 

■ arcsm — - ms m tt i (45j 



3 (1 - s)4 12/ 9(1 - s) 9 9 



.9) 87r2(4-s) 8 (5 -2s) 16 , 

)I + 97r^ + - ..j 4 + 3.-s^)arcsm 

32 (4 





+ — —-1 ^^^^^^ I ~o~ 1 ~^ cTTi ^l^'S (46) 



Fig. O shows R{s) defined in (jS]), where we set /i = 5 GeV and used y/z = mc^poic/Tnh poic = 
0.29. We used the pole-mass for the 6-quark and the MS-mass for the top-quark and set 



^^ft.poic = 4.8 GeV and mJjrit) = 163 GeV 5l|- We neglected the finite bremsstrahlung 



corrections calculated in ISj. From Fig. [3] we conclude that for /i = the contributions 
of the form factors f[^2^^ lead to corrections of the order 10% — 15% at the level of the 
perturbative part of the normalized spectrum R{s). Integrating R{s) over the high s 
region, we define 

i?high= / dsR{s). (47) 



Jo. 6 

Fig. m shows the dependence of the perturbative part of -Rhigh on the renormalization scale. 
We obtain 

i?high,pert = (0.43 ± O.Ol(^)) X 10-^ (48) 

where we determined the error by varying fi between 2 GeV and 10 GeV. The corrections 
due to the form factors -Fi lead to a decrease of the scale dependence to 2%. 

We should mention at this point that a normalization different from the one in (jH]) has 



been proposed in 2J]: By normalizing the B —>■ Xsi~^i decay rate to the semileptonic B 
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0.55 



0.45 




0.35 



FIG. 4: Perturbative part of -Rhigh as function of the renormalization scale fi at NNLL. The soHd 
line represents the NNLL result, whereas in the dotted line the order as corrections to the matrix 
elements associated with Oi 2 are switched off. See text for details. 



XuC'Ue decay rate with the same cut on g^, the large theoretical uncertainties due to power 
corrections can be significantly reduced. It was shown explicitly in a recent phenomenological 
update 261] that the uncertainties from the poorly known 0(1/ ml) power corrections are 
then under control. 



V. CONCLUSIONS 

We calculated for the first time the NNLL virtual QCD corrections of the matrix elements 
of Oi and O2 in the high region as analytic functions of and rric. While keeping the 
full analytic dependence on g^, we evaluated the matrix elements as an expansion in z 
up to the 10th power, which is numerically stable for s > 0.6. Making extensive use of 
differential equation techniques, we fully automatized the reduction of the higher order 
expansion coefficients to the leading and first subleading power, which were obtained via 
the method of regions. 

Comparing our results for these matrix elements with those of a previous work where the 
master integrals were calculated numerically we obtain an agreement up to 1%. Likewise 
in coincidence with [ij], we find that the corrections calculated in the present paper lead to 
a decrease of the perturbative part of the g^-spectrum by 10% — 15% relative to a NNLL 
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result where these contributions are not taken into account and reduce the renormahzation 
scale uncertainty to 2%. 

We provide the rather lengthy results of our calculation in electronic form as Mathematica 
files and for numerical purposes also as C++ files. 
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APPENDIX A: COMMON MASTER INTEGRALS 



All integrals are evaluated in = 4 — 2e dimensions. In the following notation we suppress 
the positive imaginary part +iO of the denominators. The integration measure is defined as 




1. One-loop integrals 

a. 2-point integral with two massive lines 




[dk] 



1 



(A2) 



{k'^ — m'^){{k + qy — m?) 
The double line denotes the massive propagator. We evaluate ( 1A2I) in the two regions 
< Am? and q^ > 4m^. In the latter one we need the integral in an expansion in m?/q^. 
Using Mellin-Barnes representation 28|, |52| it is easily seen that we can cast (lA2p into the 
following form: 



(4vr) 



q2 



1 

2Td 



dt 



m 
q2 



T{-t)T{t + e)T^{l-t-e) 
r(2 -2t- 2e) 



(A3) 



where the integration contour over t has to be chosen such that — e < 3?(t) < 0. The poles 
on the right hand side of the contour are located at t = n and t = n + 1 — e where n G Nq. 
By closing the integration contour to the right we obtain the power expansion in m? /q^. 
Now let us consider the region q^ < Aim?. Up to order e the integral reads: 

,2 




{Any 



X 



A — X \/x 
— - — arcsin h 



X 



4 + 



TT 

12 



A — X y X 

arcsin — 

X 2 



-4 + 21n(4-x)) + 2 



A- X 



Clo 2 arcsin 



X 




(A4) 



where we defined x = q^ /m? and Cl2(0) = 53Li2(e*'^ 



22 



b. 3-point integral with one massive line 




[dk] 



where 



= rn^, = xm^ and x < 1. 



The integral is evaluated to 




/iV^V 1 r(e) 



X 



(47r)2 \ m? J m? 1(2 - 2e) 
[r(l - 2e) 2Fi(l, 1; 2 - 2e; 1 - x) - x-'e'^'V\l - e) 2Fi(l, 1 - e; 2 



with 2F1 given by ( fTTj) . 



c. 3-point integral with two massive lines 



[dk] 



k'^{{k + p)"^ — m?'){{k + p — qY — rn?) 



where 



p^ = rn^, q^ = xm^, (p — qY = (f < 4m^. 
The expansion in e of flASp reads 



I ( fJ'^\^ 1 



(47r)2 \m? J m?{l—x) 



TT 2 V a; 
ho arcsm h 



e ( 12 In (1 — x) arcsin^ — h 4CI2 I 6 arcsin — h tt | arcsin ^ 



1 / vx\2/ Vx 

-TT^ In (1 — x) + 6CI3 I 2 arcsin — h tt j + ^ CI3 I 6 arcsin — h tt 



where = 3fJLi3(e*' 
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2. Two-loop integrals 



a. Two massive lines 

We need the following three sunrise diagrams in an expansion in m? /q^. So as above we 
jive the Mellin-Barnes representation, from where the expansion can be easily derived. 



[dk][dl] ^ 




r(-t)r(t - 1 + 2e)r^(l - e - t)r(2 -2e-t) 

r(2-2e-2t)r(3-3e-t) ' ^ ' 

The residues we have to take into account are located at n, n + 1 — e and n + 2 — 2e with 

n G Nn. 




2^ 



, ■ , , r(-e)— / dt [ ^ \ e*"^'^+*'x 

T(-t)Tit + 2e)T\l-e-t)Ti2-2e-t) 

r(2 - 2e - 2i)r(2 - 3e - 1) ' ^ ' 

The dotted line denotes a propagator that has to be taken squared. The residues are located 

at n, n + 1 — e, n + 2 — 2e, n G Nq. 



1 




2^ 



(fc + g)2(/2-m2)2((A; + /)2_rn 

1 ^/^Vpn_..^ r..«'.^^(2e..). 



r(-i)r(i + 2e)r(-e - t)r(l - e - t)r(l - 2e - t) 



(A13) 



r(l-2e-2t)r(2-3e-t) 
with the residues located at ra, — e, n + 1 — e, n G Nq. 

b. Three massive lines 

We need the following three integrals in an expansion in m? /M"^. Therefore we give their 
Mellin-Barnes representation. 
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[dk][dl] 



(P - M2)(/2 _ m2)2((A; + /)2 _ ^2) 
,2 \ * 



(47r)4 \ J 2m 

T{-t)T{t - 1 + 2e)r2(l - e - t)r(e + t)T{2 - 2e - t) 



r(2-2e-2t)r(2-e) 
with the residues located at n, n + 1 — e, n + 2 — 2e, n G Nq. 



(A14) 




[dk][dl] 



.2 \ i 



fj^e^y 1 p /m2Y 



(47r)4 

V{-t)V{t - 1 + 2e)r2(l - e - t)r(e + t)r(3 - 4e - 2t) 
r(2-2e-2t)r(3-3e-t) ' 

with p2 = ]Vf2 and the residues located at n, n + 1 — e, 72/2 + 3/2 — 2e, G Nq. 



(A15) 




[dk][dl] 



1 



((A; + p)2 - M2)(/2 - m2)2((A: + If - m? 



1 / /iV^^ Y 1 r°° /w_ 



2 \ * 



X 



T{-t)T{t + 2e)r(-e - t)T{l - e - t)T{l + e + t)T{l - 4e - 2t) 
r(l-2e-2t)r(2-3e-t) ' 

with p2 = ]\^2 g^j^^ ^Yie residues located at n, n + 1 — e, n/2 + 1/2 — 2e, n G Nq. 



(A16) 
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